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Abstract. We present infrared spectroscopy and imaging 
of AFGL 4106. The 2.4-5 /im ISO-SWS spectrum reveals 
the presence of a cool, luminous star {TqQ ~ 3750 K) in 
addition to an almost equally luminous F star (Tgg « 
7250 K). The 5 - 195 iim SWS and LWS spectra are dom- 
inated by strong emission from circumstellar dust. We find 
that the dust consists of amorphous silicates, with a mi- 
nor but significant contribution from crystalline silicates. 
The amorphous silicates consist of Fe-rich olivines. The 
presence of amorphous pyroxenes cannot be excluded but 
if present they contain much less Fe than the amorphous 
olivines. Comparison with laboratory data shows that the 
pure Mg-end members of the crystalline olivine and py- 
roxene solid solution series are present. In addition, we 
find strong evidence for simple oxides (FeO and AI2O3) 
as well as crystalline H2O ice. Several narrow emission 
features remain unidentified. 

Modelling of the dust emission using a dust radiation 
transfer code shows that large grains (« 1/im) must be 
present and that the abundance of the crystalline silicates 
is between 7 and 15% of the total dust mass, depend- 
ing on the assumed enstatite to forsterite ratio, which 
is estimated to be between 1 and 3. The amorphous 
and crystalline dust components in the shell do not have 
the same temperature, implying that the different dust 
species are not thermally coupled. We find a dust mass of 
« 3.9 • 10~^ Mq expelled over a period of ~ 4 • 10^ years 
for a distance of 3.3 kpc. The F-star in the AFGL 4106 bi- 
nary is likely a post-red-supergiant in transition to a blue 
supergiant or WR phase. 
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1. Introduction 

Late stages of evolution of both low and high-mass stars 
are characterized by high mass-loss in a slow and dusty 
outflow. Low- and intermediate-mass stars (up to 8 M©) 
will lose a significant amount of their initial mass on the 
Asymptotic Giant Branch (AGB). More massive stars will 
lose their hydrogen-rich envelopes during a red super- 
giant (RSG) phase. The mass-loss rates can be as high as 
10"'' M0/yr for the low-mass stars and up to 10~^ Mo/yr 
for the high-mass stars. Low and intermediate mass ob- 
jects evolve to the blue part of the HR diagram after the 
entire H-rich envelope has been exhausted, and a plane- 
tary nebula (FN) is formed. High mass stars either ex- 
plode as a supernova while they are red supergiants, or 
they evolve to the blue part of the HR diagram to become 
a population I Wolf-Rayet (WR) star before ending their 
live. 

The evolution of transition objects, i.e. stars that are 
in a rapid phase of evolution between the red and blue 
part of the HR diagram is not very well known. It is clear 
however, that mass loss plays an important role in de- 
termining evolutionary timescales. Remnants of previous 
mass-loss phases are still in the vicinity of the star in the 
form of detached dust shells. Infrared observations clearly 
reveal the existence of these dust shells and the study of 
the Circum Stellar Environment (CSE) provides essential 
information about mass-loss history and consequents of 
evolution. 

The infrared-bright 
object AFGL 4106 (= IRAS 10215-5916) was classified 
as a transition object by Hrivnak et al. (1989, hereafter 
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HKV) based on its IRAS colours. Some confusion exists in 
the literature about its identification. HKV quote Bidel- 
man, who claim that IRAS 10215-5916 is equal to CoD - 
58°3221. Apparently, this has sometimes been mistaken in 
the literature as CPD -58°3221 and therefore HD 303822 
(e.g. in the Tycho catalogue). The correct identifications 
should be HD 302821 and CPD -58°2154. 

The central star was classified as a G2 supergiant by 
Garcia-Lario et al. (1994). The combination of a cool cen- 
tral star with low surface gravity and a relatively warm 
(i.e. young) detached dust shell is typical for transition 
objects. Garcfa-Lario et al. (1994) found double-peaked 
[N II] lines in the optical spectrum, which they thought to 
originate from an expanding envelope. The presence of the 
[N II] is surprising given the spectral type of this object. 
From the separation between the two peaks they found an 
expansion velocity of Voxp = 17±2 km s~^, which is quite 
reasonable for an AGB-wind. The expansion velocity of 
the detached shell was also measured in the CO (J= 1 — 0) 
rotational line and was found to be about 40 km s~^ (Jos- 
selin et al. 1998). This is high for a low mass AGB star 
and suggests a massive progenitor. The difference with the 
[N II] line is remarkable and is probably related to the ex- 
tended nature of the [N II] emitting region (see van Loon 
et al., in prep.). 

The distance and luminosity of AFGL 4106 are un- 
certain, which makes it hard to decide directly whether 
the progenitor was an AGB star or a RSG. The luminos- 
ity of K, 10^ X dP (with d in kpc), suggests that 
it is a luminous object, which until 2 kpc is a post-AGB 
star, however if this object is located at a larger distance 
a post-RSG nature cannot be excluded. We will present 
evidence that AFGL 4106 is likely a post-RSG. The op- 
tical to sub-millimetre energy distribution of AFGL 4106 
was analysed in detail by HKV and Giirtler et al. (1996, 
hereafter GKH). Both studies find a very high mass loss 
which recently stopped, however with rather different den- 
sity distributions in the shell. 

AFGL 4106 was selected for observations with the In- 
frared Space Observatory (ISO) in the guaranteed time 
of the Short Wavelength Spectrometer (SWS) consortium 
because of its classification as a transition object, and 
because of its very high infrared flux levels. The high- 
quality ISO data allow a detailed determination of the 
dust composition and its distribution in the circumstellar 
shell, which may yield new insight in the origin and evolu- 
tion of this remarkable object. Waters et al. (1996) studied 
the 30-45 /xm SWS spectrum of AFGL 4106 and discov- 
ered weak narrow emission bands that can be attributed 
to crystalline silicates. Here we present a full inventory 
of the solid state emission bands and their identification. 
For the first time a full radiative transfer model of a CSE 
including crystalline silicates is presented. 

This paper is organised as follows; In Sect. 2 we re- 
port on the observations and we discuss the data reduc- 
tion procedures. In Sect. 3 we present the stellar near in- 
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frared spectrum (2.4 - 5 /im), which shows evidence for 
the binary nature of this object. In Sect. 4 the infrared ex- 
cess at 5 ~ 195 /im, caused by the dust shell, is analyzed. 
We identify the dust composition using the different solid 
state features detected in this part of the spectrum. In 
Sect. 5 we present a radiative transfer fit to the spectrum. 
In Sect. 6 we discuss the results of this paper and in Sect. 7 
we summarize our main conclusions. 

2. The observations and data reduction 

2.1. ISO observations 

Spectra of AFGL 4106 were obtained using the Infrared 
Space Observatory (ISO, Kessler et al. 1996). Three spec- 
tra were taken with the Short Wavelength Spectrometer 
(SWS, de Graauw et al. 1996) and one with the Long 
Wavelength Spectrometer (LWS, Clegg et al. 1996, Swin- 
yard et al. 1996). An overview of the ISO observing log 
can be found in Table [l[ Part of spectrum II has already 
been shown in Waters et al. (1996). Because of its better 
signal to noise, we will use spectrum IV for our discussion 
in this paper, but information from the other spectra has 
been used in the reduction of spectrum IV. 

2.1.1. SWS 

The SWS spectra (2.36 ~ 45.3 nm) consist of 12 sub- 
spectra with slightly overlapping wavelength coverage and 
each observed by 12 detectors. Each spectrum is scanned 
twice, the so called up- and down-scan (de Graauw et al. 
1996). The spectra were reduced using the SWS off-line 
analysis pipeline, version 6.0. For a description of flux and 
wavelength calibration procedures, we refer to Schaeidt et 
al. (1996) and Valentijn et al. (1996). In the reduction 
process the calibration files available in April 1998 were 
used. Judging from the overlap regions of the different 
sub-bands, further improvement is still possible, but will 
not influence the overall conclusions of our results. 

The main fringes in the 12.0 - 29.5 /im part of the 
spectrum were removed using the InterActive (lA) data 
reduction package routine fringes. Major irregularities due 
to glitches and large drops or jumps in each detector were 
removed by hand. A comparison between the different de- 
tectors scanning the same wavelength region was used to 
determine the location of such jumps. Detectors with a 
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poor performance at the time of observations were re- 
moved. 

Several procedures were performed on each sub-band 
separately to combine the 12 detectors. A sigma clipping 
procedure was used to remove deviating points. The flux 
level of the 12 detectors was brought to the average level of 
the 12 detectors. Finally each sub-spectrum was rebinned 
to a resolution of A/AA = 750. 

The individual sub-spectra, when combined into a sin- 
gle spectrum, can show jumps in flux levels at band edges 
due to imperfect flux calibration or dark current subtrac- 
tion. Based on the observed jumps, downwards with in- 
creasing aperture size, we assume that the flux differences 
are not due to the significant difference in beam size. In- 
deed, our ground-based 10 //m image (see Sect. 2^) shows 
an angular size well within the SWS apertures. We have 
adjusted the different subbands, according to the expected 
source of discrepancy, to form a continuous spectrum. 
Note that all the shifts are well within the photometric 
absolute calibration uncertainties (Schaeidt et al. 1996). 
The 29.5 - 31 /im section is severely affected by memory 
effects of the detectors. Therefore we have removed this 
part from the final spectrum. 

In Fig. |l| we show the final SWS spectrum and we 
compare this to the original (without band matching) data 
of spectrum IV. A comparison of the final SWS spectrum 
with the original data of spectrum II gives a match of 
similar quality. 



2.1.2. LWS 

Spectrum III is an LWSOl full spectral scan (43 - 195 /im). 
The spectrum was sampled with one point per spectral 
resolution element using fast scanning. An off-set position 
spectrum was also taken with the same instrumental pa- 
rameters, in order to correct for background contributions. 
The spectrum was reduced using the LWS off-line analysis 
pipeline, version 6.0. See Table |l| for an overview of the 
observing mode and date of observation. 

The Auto Analysis results of both source and back- 
ground spectrum were then processed using the ISO Spec- 
tral Analysis Package (ISAP). First a median clipped av- 
erage was calculated of all the spectral scans for one de- 
tector. These averaged scans were then fitted together us- 
ing the scan on the LWl detector as reference. The re- 
sulting spectrum was rebinned to a standard wavelength 
scale with a resolution of 250 using a flux conserving re- 
binning routine. Finally, the on-source and off-source spec- 
trum were subtracted. The final SWS + LWS spectrum is 
presented in Fig. |[ 

2.1.3. Comparison with IRAS data 

In Fig. ||we show the complete spectrum together with the 
IRAS broad band fluxes. We compared the resulting 2.4 - 
195 fim spectrum with IRAS data by convolving the ISO 
spectrum with the IRAS photometric band transmission 
curves (IRAS explanatory supplement 1988), the resulting 
fluxes are given in Table I 
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Fig. 1. The final SWS spectrum (thick line), with the dif- 
ferent sub-bands attached to each other, compared with 
the original outcome of the data reduction (dots). 



Band IRAS-flux Folded ISO-flux 



{/iin) (Jy) (Jy) 

12 200.8 169 

25 1755 1883 

60 851.8 910 

100 181.1 283 



Table 2. The comparison with the IRAS broadband flux 
density 



The differences found for the 25 and 60 /im flux den- 
sities are well within the error that IRAS, SWS and LWS 
quote for their flux uncertainty. This also holds for the 12 
/im flux density, however one should realize that discrep- 
ancies in this wavelength regime might be related to the 
memory effects found in band 2C both for this source (as 
was found as a difference in the up- and down-scan) as 
well as in the calibration sources. 

The 100 /im IRAS flux density (181.1 Jy) is consid- 
erably lower than the flux density in the LWS spectrum 
(283 Jy) and is outside the errors quoted for IRAS and 
LWS. We have inspected the raw IRAS data covering the 
position of AFGL 4106 for bad scans or steep background 
gradients, which may cause the difference between both 
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observations. We found no bad scans but there is a con- 
spicuous gradient in the galactic background radiation. 
However, the IRAS 100 /im flux is corrected for this back- 
ground. Although the contribution of the background to 
the LWS continuum flux was found to be minor, the pres- 
ence of the [C II] line at 158 /xm, likely to be of interstel- 
lar origin, in the continuum subtracted spectrum suggests 
that the background has not been subtracted correctly. At 
present, the difference between IRAS 100 iim and LWS is 
not understood. In Fig. ||, we also show the comparison 
of the IRAS Low Resolution Spectrograph (LRS) and the 
final spectrum. The agreement of the shape and flux levels 
of the LRS and the SWS spectra is satisfactory. 




5 10 50 100 
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Fig. 2. The final SWS and LWS spectrum of AFGL 4106 
(solid line) compared with the IRAS-LRS (dotted line) 
spectra and the IRAS broad band fluxes (triangles) 



2.2. 10 /im imaging at the ESO 3.6m telescope 

On 23 February 1994 we used the mid-infrared camera 
TIMMI (see Kaufl et al. (1994) for a description of the in- 
strument) at the ESO 3.6 m telescope on La Silla, Chile, 
to image AFGL 4106. We used a broad band N filter, 
centered at 10.1 /im. Standard infrared techniques were 
used, i.e. chopping and nodding with a throw of 23.5" , 
somewhat larger than the 21.55" field of view. We inte- 
grated for 31 minutes on the source. We chose a small 
pixel scale of 0.336" to avoid undersampling of the Point 
Spread Function (PSF). We observed the standard stars 
a Hya, e Mus and a Car, to allow for absolute flux calibra- 



tion and to derive the PSF (FWHM w 1.2") of the tele- 
scope plus instrument. We reached the diffraction limit of 
0.6" by applying a deconvolution algorithm (Lucy 1974) 
within the MIDAS reduction software package. This de- 
convolution method has as drawback that it focuses the 
noise in the background, which is therefore quite spiked. 
The deconvolved image is shown in Fig. pi The integrated 
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Fig. 3. A logarithmic scale plot of the N-band image of 
AFGL 4106. The black line contours start dX <j — 1 and 
are steps of factors 3, the white line start at ct = 300 and 
increase in steps of 100 up to 1000. North is up and East 
is to the left. 

instrumental 

N-band magnitude (not corrected for color differences 
etc.) is iV -1.19 ± 0.05. This corresponds to 120 Jy. 
The convolution of the N-band sensitivity with the ISO 
spectrum gives a flux of 100 Jy. The difference between 
both flux levels is probably caused by detector memory 
effects in the ISO band 2C wavelength region, and by the 
slightly extended nature of the nebula. The brightest, un- 
resolved emission peak in the image is likely to be centered 
on the central object. The nebula, especially the outer 
parts, seems to have an oval to box- like shape. The inner 
part of the emission complex is characterized by structure 
on (sub-)arcsecond scale in the form of local enhancements 
of 10 /im emission. Whether this also reflects a clumpy 
density structure is not certain. The North- West part of 
the nebula is fainter compared to the rest of the nebula 
at 10 /im, and seems slightly more extended. The whole 
nebula and especially the inner part has a North- West - 
South-East symmetry axis. The unresolved central peak 
contains roughly 13% of the total flux, which is about 15 
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Jy. The total emission extends out to « 5" from the cen- 
ter. 

An arc of Ha emission, extending from roughly North- 
East to South- West clockwise, has been found (van Loon 
et al., in prep.), and is located just outside the edge of the 
extended dust emission in the N-band. This arc coincides 
with the low surface brightness emission region in the N- 
band image. 



3. The 2.4 

revealed 



5 fim part of the spectrum; A binary 



The spectral energy distribution (SED) can naturally be 
divided into two parts. For wavelengths shorter than 5 
/zm the photosphere radiation dominates the spectrum. 
Above 5 /im the thermal radiation of the dust dominates 
the spectrum. This IR-excess will be discussed in Sect. ^. 
Here, we will present evidence that the central object is 
a binary system, consisting of an M and A-F type star of 
almost equal luminosity. 

The 2.4 — 5/Ltm part of the spectrum is shown in Fig. ^. 
It is dominated by strong absorption from several band 
heads of ro-vibrational molecular bands, such as the CO 
first overtone and fundamental, the SiO first overtone, 
H2O (vi and 1^3 ) and OH fundamental. The strength and 
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Fig. 4. Ro-vibrational band heads of different molecu- 
lar species, proving the existence of a cool companion to 
HD 302821 

excitation levels of these bands are not compatible with 
those expected for a G-type supergiant spectrum, but sug- 
gest a cooler star. From the shape and the level of the first 



overtone ro-vibrational ^^CO bands we derived a tempera- 
ture of « 4000 K (see e.g. a-Boo in Decin et al. 1997). This 
high temperature together with the required column den- 
sity cannot be caused by circumstellar gas in a detached 
(dusty) envelope. Also the presence of the high excita- 
tion lines of gas-phase SiO, which cannot be in a dusty 
envelope, suggests that we are dealing with a cool star. 
Meanwhile van Winckel et al. (in prep.) took a high res- 
olution spectra of AFGL 4106 in the optical and around 
1/im and found a warm star of T w 7000 — 7500 K, domi- 
nating in the optical and clear evidence of both the warm 
and the cool star around 1 /zm, where both stars appear 
to be equally luminous. 

In order to verify the binary hypothesis, we have com- 
bined the optical and near-IR photometry of AFGL 4106 
(HKV, Garci'a-Lario et al. 1994, 1997) together with our 
ISO data. We used Kurucz (1991) model atmospheres to 
fit the broad-band photometry up to 2.2 /im. This ensures 
that the fit is not contaminated by dust emission. For our 
fitting procedure we assumed that the mean interstellar 
extinction law can be used to describe both the circum- 
stellar and the interstellar dust extinction. We will come 



back to this point in Sect. 5.2. We were not able to fit 
the data with a single star, but found a good fit if we as- 
sumed a binary. A best fit was found with temperatures of 
7250±250K and 3750±250K for the hot and cool compo- 
nent respectively, while the luminosity ratio (Lwarm/Lcooi) 
= 1.8. For both components we independently found the 
same reddening of E{B — V) ~ 1.2 ± 0.1 mag. The fit 
together with the observations is shown in Fig. 

The relatively small luminosity difference between 
both stars suggests that both components are evolved, i.e. 
the M star is a giant or supergiant. A more detailed com- 
parison of the ISO 2.4-5 fim spectrum with that of other 
M giants and supergiants indicates a supergiant nature of 
the spectrum (Van Winckel et al., in prep). 

The flux from the central part of the ground-based 10 
/im image is about 15 Jy. This is much larger than the 
flux expected from the stellar atmospheres of both stars 
in the system (3 Jy) . This suggests that there is additional 
circumstellar material very close to the binary, probably 
hot dust. It is not clear whether this material was ejected 
at the same time as the dust and gas in the detached shell, 
or whether this represents a later phase of mass loss. The 
additional circumstellar material could also originate from 
the M-star in the system. Note that the flux of the second 
brightest blob, SW of the central peak, is still about 7 Jy, 
which is also enough to contain one or even both stars. 
However one should keep in mind that both stars show the 
same reddening, suggesting that they are located behind 
the same dusty material and that more clumps (e.g. SE, 
E and NNE of the central peak) are present. 



6 



F.J. Molster et al. 



■10.5 



DO 

o 



11.5 




12.5 



Fig. 5. The reddened {E{B — V) — 1.2) composite spec- 
trum (thick hne) of the hot star (dotted line) and the cool 
star (dashed hne) compared with the broadband photom- 
etry (triangles) and the ISO-data (thin line). The error 
in the photometry is much smaller than the size of the 
triangles. 

3.1. Interstellar and circumstellar extinction 

The large reddening of AFGL 4106 is illustrated by the 
prominent diffuse interstellar bands (DIB's) in the optical 
spectrum of AFGL 4106. In Fig. || we display part of the 
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Fig. 6. High-resolution optical spectrum of AFGL 4106. 
The prominent interstellar DIB's are indicated together 
with the main photospheric atomic absorption lines. The 
spectrum has been corrected for the presence of telluric 
lines. 

optical spectrum taken by van Winckel et al. (in prep.) 
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Table 3. Equivalent width of a selection of DIB's in the 
optical spectrum of AFGL 4106 compared to the strength 
of same DIB's in BD-|-63°1964 (Ehrenfreund et al. 1997) 
and the DIB's- reddening relation of Jenniskens & Desert 



corrected for telluric lines. The broad band at 6284 A, 
which is a blend of different components, is strong, but 
also the narrow DIB's at 6376 and 6379 A can easily be 
recognized. Since the carriers of the DIB's are associated 
with the carbon-rich component of the interstellar medium 
(ISM) and the circumstellar material around AFGL 4106 
is oxygen-rich, the DIB's are likely to be purely from in- 
terstellar origin and can be used to estimate the inter- 
stellar component of the total reddening of AFGL 4106. 
The strength of the DIB's scales only roughly with the 
reddening {E{B — V)) (see Herbig 1995), therefore our 
quantisation of the interstellar reddening is only indica- 
tive. 

We compared the strength of the DIB's in AFGL 4106 
with the DIB strength in BD+63°1964 {E{B ~ V) ^ 
1.01 mag.. Spectral type BOII) (Ehrenfreund et al. 1997) 
and derived a second estimate of the reddening by us- 
ing the equivalent width to reddening relation derived by 
Jenniskens & Desert (1994) from a detailed analysis of 
4 reddened early-type stars. We limited ourselves to nar- 
row DIB's without photospheric contamination since the 
broad DIB's in our spectrum are composite and the equiv- 
alent width determination most inaccurate. The equiva- 
lent width measurements are given in Table | One can see 
that the DIB's in AFGL 4106 and BD-f63°1964 are very 
similar. We estimate the reddening using a linear least 
mean square analysis and found E{B — V) = 1.07 mag. 
by using the mean relation of Jenninskens & Desert and a 
similar value E{B-V) = 0.94 mag. by using BD-h63°1964 
as a standard. We conclude that the interstellar reddening 
towards AFGL 4106 is E{B - V) = 1.0 ± 0.2 mag. This 
is in very good agreement with the value derived by van 
Loon et al. (in prep). 

The case of AFGL 4106 offers a unique opportunity 
to measure the internal extinction within the cicrumstel- 
lar dust shell. This is because the star and the dust shell 
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are surrounded by a shell of ionized gas, detected in Ha 
and the [N II] line. These emission lines were reported 
by Garcia Lario et al. (1994), and subsequent imaging 
of the spatial extent of the emission lines by Van Loon 
et al. (in prep.) showed that the lines originate from a 
roughly spherical shell of gas on the periphery of the dust 
shell. The [N II] lines were found to be split, with the red- 
shifted component systematically weaker than the blue 
peak. Van Loon et al. (in prep.) interpret this difference 
in intensity as a result of internal extinction in the dust 
shell, and assuming a reddening law which is the same ' 
for the interstellar as well as for the circumstellar dust, 
they found a circumstellar reddening towards the cen- 
tral star of E{B - V) = 0.22 ± 0.05 mag. We conclude 
that sum of the DIB (interstellar) extinction of about 
E{B — V) = 1 mag. and the [N II] (circumstellar) extinc- 
tion of about E{B — V) = 0.2 mag. agrees well with the de- 
rived total line-of-sight extinction of E{B — V^) = 1.2 mag. 
based on the SED. 

From the total spectrum one can derive the total flux, 
infrared excess and therefore another estimate of the inter- 
stellar and circumstellar absorption towards AFGL 4106. 
We define the infrared excess as the integrated flux under 
the energy distribution corrected for the stellar contribu- 
tion. The distance to AFGL 4106 is not known, so we can 
only determine a distance dependent luminosity, which is 
roughly 10.8 x 10'^ rf^L©, with d the distance in kpc. The 
IR-excess contributes about 75% (8.0 x lO'^d^L©) to the 
total luminosity. Which leaves 2.8 x 10^ (PLq for the lu- 
minosity directly from the reddened photospheres of both 
stars. Assuming an E{B — V) of 1.2 mag. for both stars, 
taken from the fit to the broad-band photometry, results 
in a dereddened luminosity of ~ 7.4 x lO'^d^L© for the 
hot component and w 4 x 10'^ (PLq for the cool one. This 
implies a total "loss" of stellar energy by absorption and 
scattering of « 8.6 x lO'^d^Lo/s. This would imply that, 
assuming a spherical dust shell and assuming that all ab- 
sorbed energy is re-radiated in the infrared, more than 
90% of the extinction towards AFGL 4106 is circumstel- 
lar, and thus the contribution of the interstellar extinction 
would be negligible. This is not in agreement with the esti- 
mate of the interstellar extinction derived from the depth 
of the DIBs (see above). In Sect. 5.2 we will demonstrate 
that this discrepancy is resolved by the non standard ex- 
tinction properties of the CSE. 

4. The 5 195 /im part of the spectrum; The 
dust shell 

At wavelengths longwards of about 5 /im thermal emission 
from the circumstellar dust starts to dominate the spec- 
trum (see Fig. Between 5 and 7.8 fim the spectrum is 
flat at a level of approximately 15 Jy. Around 7.8 fim a 
steep rise is found with a flattening around 12 ^m. This 
prominent structure is caused by an increasing continuum 
together with the 10 /J,m amorphous silicate feature (Si-0 



2000 




50 
A [/u.m] 

Fig. 7. The SWS and LWS spectrum of AFGL 4106 with 
the spline fit continuum. The result of the subtraction of 
the spline fit continuum (dashed line in A) is shown in 
graph B. The result of the removal of the broad silicate 
features (dashed line in B) is given in graph C. 



stretching mode). At 13 /im the flux increases again, with 
at around 18 /im the second peak of the amorphous silicate 
feature (0-Si-O bending mode). The peak of the infrared 
excess is found around 25 /im. Superimposed on this con- 
tinuum with the broad features are narrow features, es- 
pecially in the SWS part of the spectrum. In most cases 
they can be attributed to crystalline silicates. 

In contrast with the SWS spectrum with its sharp 
spectral features, the LWS spectrum of AFGL 4106 is very 
smooth. Emission lines are present in the on-source spec- 
trum at 88.356 ([O III] line) and 157.741 fim ([C II] line). 
However, since they are also present in the background 
spectrum, they cancel out almost completely in the final 
spectrum. Although the warm component in the binary is 
not hot enough to produce the [O III]-line, the presence 
of the [O III]-line may not be surprising given the line of 
sight (Carina) and the projected vicinity of several very 
luminous hot stars, which emit enough ionizing radiation 
to produce the [O III]-line. We find no evidence for emis- 
sion lines in the SWS and LWS infrared spectra from the 
ionized nebula surrounding AFGL 4106 (detected at opti- 
cal wavelengths), which should be within the ISO beam. 
Clearly the dust continuum swamps any line emission from 
that region. 
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4.I. The continuum subtraction of the SWS and LWS 
spectra 

In order to show the presence of the spectral features more 
clearly, we have subtracted a continuum by fitting a spline 
curve to the spectrum. Note that the spline fit continuum 
has no physical meaning and is only used to enhance the 
sharp solid state features on top of the spectrum. The con- 
tinuum subtracted spectrum is shown in Fig. panel B. 
The amorphous silicate features dominate the continuum 
subtracted spectrum below 20 /im. In order to show the 
narrow bands more clearly, we have created a second con- 
tinuum subtracted spectrum, where we have excluded the 
broad amorphous silicate bands (see the dashed line in 
Fig. panel B and the result in panel C). We find a wealth 
of solid state structures, which we show in more detail in 
Fig. ||. Some residual fringing is still visible at 18-20 /xm. 

Since the solid state bands have a wide range in width 
and strength, and often blend, we decided to measure the 
band strength with respect to a local continuum. We used 
a second order polynomial fit to define a local continuum 
in the final spectrum (i.e. not the continuum subtracted 
spectrum shown in Fig. ^and^), and fitted Gaussians to 
the emission bands with respect to this local continuum. 
The very broad amorphous silicate bands could not be 
fitted with a single Gaussian, and we give the FWHM 
directly measured from the spectrum. The results of the 
Gauss-fits can be found in Table ^. Due to our continuum 
subtraction method, it is likely that very broad features 
have been included in the continuum and therefore are not 
visible in the continuum subtracted spectrum. 

In the following sub-sections, we briefly describe the 
continuum subtracted spectrum. 

4.1.1. The 10-20 /^m region 

The Si-0 stretching and 0-Si-O bending mode vibrations 
of the amorphous silicates dominate the 10 - 20 /im part 
of the spectrum. We note that the peak wavelength of the 
Si-0 stretch is significantly red-shifted compared to the 
bump found in most 0-rich (post-)AGB stars (Monnier 
et al. 1998) and the ISM. However, Little-Marenin and 
Little (1988, 1990) find in LRS spectra an emission feature 
which peaks at a similar position in most S-, several MS- 
and some (5%) M-type stars. On top of the amorphous 
silicates are narrow emission features. Careful inspection 
of the raw data suggests that there are three instrumental 
artifacts at 11.07 and probably 13.6 and 14.2 /im, but 
the other features at 16.1, 16.8, 18.1 and 19.2 /im are 
of circumstellar dust origin. The identification of most of 
these narrow peaks can be found in Table H and Sect. 4.2. 



tures with a range of FWHM; we find narrow features 
(FWHM/A < 0.03) at 20.6, 21.5, 22.9 and 24.1 /im, and 
broader features (FWHM/A > 0.03) at 23.6, 26.2, and 
27.8 /im. The sharp descent at the long wavelength side 
of the 21.5 /im feature and rise at the short wavelength 
side of the 22.9 /im feature are probably caused by an 
imperfect responsivity calibration. 

4.1.3. The 29.5-45 /im region 



As stated in Sect. 2.1.1 we excluded 29.5 to 31 /im from 



4.1.2. The 20-29.5 /im region 

The spectral region between 20 and 29.5 /im is shown in 
Fig. [s] and is characterized by a wealth of emission fea- 



the final spectrum, but we can say that we did not detect 
narrow emission bands in this wavelength range. The spec- 
tral region from 29.5 to 45 /im is dominated by a broad 
structure from about 31 to 44 /im, on top of which emis- 
sion features are present, see also Fig. |. We find broad 
features at 33.6, 40.4 and 43.1 /im and narrow features at 
32.8, 34.0, 35.8, 36.5 and 41.1 /im. Waters et al. (1996) find 
a weak plateau from 32 to 37 /im in the spectrum of sev- 
eral oxygen rich dust shells. This plateau is also visible in 
AFGL 4106, although it is weak and blends with the broad 
unidentified structure from 31 to 44 /im. There is also evi- 
dence for sub-structure in the 32-33 /im region, suggesting 
another weak band that contributes. This weak band at 
32.0 /im is also seen in the SWS spectrum of NGC 6302 
(Lim et al., in prep.) and in that of HD44179 (the "Red 
Rectangle", Waters et al. 1998). 

4.1.4. The 45-195 /im region, LWS 

At longer wavelengths, the spectrum is smooth. We find 
a narrow band at 47.8 /im, and a broad emission centered 
at 61 /im (probably due to H20-ice). Also, the [C II]-line 
at 157 /im is present; this is likely due to imperfect can- 
cellation of the on-source and off-source caused by the 
gradient found in the background (see also the remarks 
in Sect. 2.1.3 ). For the moment, we ignore the [C II]-line 
emission. 



4. 2. Identifications 

4.2.1. Amorphous silicates 

The broad emission features around 10 and 18 /im are at- 
tributed to amorphous silicates. We have derived a rough 
estimate of the temperature of 140 K for the amorphous 
dust component using the amorphous olivine optical con- 
stants published by Dorschner et al. (1995) (see also 
Fig. ^) . Such a temperature estimate should be used with 
caution, since the features are on top of a steep contin- 
uum and contribution from other species (e.g. FeO) can 
affect the 10 over 18 /im band strength ratio. A similar 
temperature was found by fitting the shape of the dust 
spectrum with a dust emissivity Q(A) oc A^^. This simple 
procedure even fits the long wavelength part very well but 
is significantly too low from 15 to 50 /im. 
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Table 4. The characteristics of the dust features in the spectrum of AFGL 4106 



4.2.2. Crystalline olivines and pyroxenes 

We use the laboratory transmission spectra of crystalline 
silicates as measured by Jagcr ct al. (1998) and Koike & 
Shibai (1998) to identify the rich solid state spectrum. The 
identifications are based on the match between the labo- 
ratory peak positions and observed peak positions for the 
minerals of the olivine (Mg2a;Fe(2_2x)Si04) and pyroxene 
(Mg2:Fe(i_a;)Si03) family. We find a very good match be- 
tween the most (but not all) broad emission bands and the 
laboratory data for the pure Mg end member of the crys- 
talline olivine solution serie, i.e. forsterite (Mg2Si04). Sim- 
ilarly, we can identify most narrow emission bands with 
the pure Mg end member of the crystalline pyroxene so- 
lution serie, i.e. enstatite (MgSiOs). Enstatite comes in 



two crystal structures, ortho- and clino-enstatite. How- 
ever, apart from the 32.8 /im peak which is only found in 
ortho-enstatite (Koike & Shibai 1998) and blends with the 
33.6 /zm forsterite peak, the main strong peaks of clino- 
and ortho-enstatite are located at about the same posi- 
tion. For this reason we were not able to distinguish be- 
tween these two forms in Table ^. 

Apart from the wavelength of the bands, it is also im- 
portant to consider the relative strength of the emission 
bands in the identification process, since it gives a rough 
estimate of the temperature. To estimate the tempera- 
ture we compared the continuum subtracted spectrum of 
AFGL 4106, with the continuum subtracted spectrum of 
forsterite (and enstatite and diopside in Fig. |l^). The con- 
tinuum subtraction has been performed in the same way 
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Fig. 8. The 10 to 44 /im region of the continuum subtracted spectrum enlarged. Note the richness of sohd state 
features. 



as it was done for AFGL 4106 (i.e. a smooth sphne fit 
curve without obvious kinks) after multiplying the emis- 
sivity of the different dust species with blackbodies of the 
given temperatures. Since the peak of the features are 
at same wavelength for the two temperatures, which had 
quite different continua, it shows that this continuum sub- 
traction does not shift the peaks. For forsterite we show in 
Fig. 1^ the best fit to the continuum subtracted spectrum of 
AFGL 4106, assuming that the forsterite grains have a sin- 



gle temperature. The best fit temperature is about 80 K. 
We stress that these temperature estimates are only ap- 
proximate; a full treatment of the radiative transfer is dis- 
cussed in Sect. |[ The temperature of the amorphous dust 
is significantly higher than that of the crystalline com- 
ponent: a temperature of 140 K for the forsterite would 
have resulted in significant emission in the 11.2 fim band, 
which is not observed. We conclude that amorphous and 
crystalline dust are not thermally coupled and that they 
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Fig. 9. The emission spectrum of a 120, 140 and 160 K grain of amorphous ohvine (Mgo.8Fei.2Si04, Dorschner et aL 
1995) compared with the spectrum of AFGL 4106 (upper panel), and the continuum subtracted emission spectrum of 
a 80 and 140 K crystaUine forsterite (Jager et al. 1998) grain compared with the continuum subtracted spectrum of 
AFGL 4106 (lower pannel). 



either have quite different optical properties, or are spa- 
tially distinct. 

Unfortunately, the situation is rather disappointing for 
the pyroxenes. The observed band strengths of AFGL 4106 
poorly match the band strengths seen in the laboratory 
data published by Jager et al. (1998). Therefore we can- 
not derive a reliable characteristic temperature for the py- 



roxenes, based on these laboratory data. We have also 
compared the band strengths of pyroxenes measured by 
Koike & Shibai (1998) to our AFGL 4106 data which 
gives a slight improvement of the match, especially in the 
32.8 /im range (see Fig. ^0|). A temperature of 80 K gives 
again a reasonable fit to the data. The difference between 
these two laboratory spectra might be caused by slight dif- 
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Fig. 10. Upper panel: The continuum subtracted emission spectrum of a 80 K crystalline ortho-enstatite grain from 
Koike & Shibai (1998) (thin solid line) and from Jagcr ct al. 1998 (dashed line) compared with the continuum 
subtracted spectrum of AFGL 4106 (thick solid line). Lower panel: The continuum subtracted emission spectrum of a 
80 K crystalline diopside grain from Koike & Shibai (1998) (thin solid line) and a diopside measured by Jager (private 
communication) compared with the continuum subtracted spectrum of AFGL 4106 (thick solid line). 



ferences in the chemical composition. The natural ortho- 
enstatite sample of Jager et al (1998) is not the pure Mg 
end member of the crystalline pyroxene solution series but 
contains 2.1 mass % FeO, while the data of Koike & Shibai 
(1998) is of a synthetic sample of enstatite with some small 
impurities from Li02, M0O3 and V2O5. Wheter this dif- 



ference at 32.8 /xm is due to the small amount of FeO, 
which has influence on the peak over continuum ratio (see 
e.g. Fig. 2 and Fig. 3 of Jager et al. (1998)), or due to the 
metaloxides found in the ortho-enstatite of Koike & Shibai 
(1998) is not known yet, but we will tentatively identify 
it with enstatite. For both the olivines and the pyroxenes 
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the observed FWHM of the bands are significantly smaller 
than those measured in laboratory spectra. It is interest- 
ing to note that both in the laboratory spectra and in our 
ISO data the olivine bands are broader than the pyroxene 
bands. 

The main component for the plateau from 32 to 37 
/im is unidentified yet. Its appearance and shape are in- ^ 
fluenced by the forsterite bands at 33.6 and 36.5 fim and'—' 
the enstatite band(s) at 35.8 (and 32.8) fim, which makes p 
it difficult to extract the unblended shape of this feature ^ 
and therefore its identification. 

4.2.3. Diopside 

A third silicate that was checked for its presence is diop- 
side (CaMgSiOa). In Fig. |l^ we show the comparison of 
diopside with the spectrum of AFGL 4106. Although the 
two datasets (Koike & Shibai, 1998; Jager, private com- 
munication) used for comparison show some differences, 
probably due to a somewhat different chemical composi- 
tion, both datasets predict prominent emmision at 21 and 
25 /Ltm, which is not observed in our spectrum. We con- 
clude that the presence of diopside cannot be confirmed 
at present. 

4.2.4. Aluminum rich dust 



Assuming that the 43.1 /im feature is a blend of different 
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In Sect. 5.2 we will demonstrate that corundum (a-Al203) 
is present in the dust shell of AFGL 4106. According to the 
condensation scheme published by Tielens (1990) corun- 
dum can transform into melilite (Ca2Al2Si07). It is not 
known by the authors if melilite will form below or above 
the glass temperature. Since we only have optical data for 
amorphous melilite (Mutschke et al. 1998) we will concen- 
trate on this. Amorphous melilite shows broad features at 
10, 15, 19 and 34 /xm. Because of the amorphous structure 
the features are very broad and difficult to identify in the 
spectrum. However the very broad (several /im) feature at 
34 fim might be a good candidate for the broad structure 
from about 31 to 44 /im found in our ISO-spectrum. 

4.2.5. Water ice 

We have identified the feature at 61 /im, hereafter referred 
to as the 60 /zm feature, with crystalline H2O ice. Omont 
et al. (1990) identified a similar feature observed with 
KAO in IRAS 09371-^1212 (Frosty Leo) with crystalline 
H2O ice, while Barlow (1998) reports the presence of 60 
/im emission in several dusty circumstellar shells. Labo- 
ratory data published by Bertie et al. (1969) show that 
crystalline H2O ice has two strong bands, one at about 43 
/im, the other at roughly 61 /tm. In the SWS spectrum 
we already found a broad and strong emission feature at 
43.1 /im, which is probably a blend of crystalline enstatite 
and crystalline H20-ice. The fact that a 60 /im feature is 
present strengthens the crystalline Il2 0-ice identification. 
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Fig. 11. The emission spectrum of a 30, 40 and 50 K crys- 
talline H20-ice grain (Bertie et al. 1969) compared with 
the continuum subtracted spectrum of AFGL 4106. 

features, including crystalline H2 0-ice, we have derived an 
upper limit of the temperature of only 40 K. This is again 
much lower than the amorphous silicate temperature but 
also significantly lower than the crystalline silicate tem- 
perature. It cannot be excluded that the 60 /tm feature 
is in fact a blend. We note that the laboratory spectra of 
clino-enstatite (Koike & Shibai 1998) show a weak 60 /tm 
bump. This would result in a higher temperature for the 
crystalline H2O ice. Still it points to another grain popu- 
lation which is not thermally coupled with the others. 

4.2.6. Iron oxides 

The identification of Mg-rich and Fe-poor crystalline sili- 
cates in non-negligible abundances suggests that some Fe 
may be present in the form of simple oxides. For spheri- 
cal FeO grains the most important infrared feature long- 
wards of 10 micron is a relatively sharp feature at 20 /im, 
which is not present in our spectrum. However FeO is very 
sensitive to shape effects (Henning & Mutschke 1997). In 
the case of a continuous distribution of ellipsoids (CDE, 
Bohren & Huffman 1983) the 20/tm peak broadens signifi- 
cantly, the FWHM becomes roughly 8 /tm, and also shifts 
to 24 /im. This feature is temperature dependent (Henning 
& Mutschke 1997) but only in absolute strength and not 
in wavelength. The broad feature is present in our total 
spectrum, but hidden in the "continuum" in the contin- 
uum subtracted spectrum. If FeO is present it is possible 
that Fe203 (hematite) and Fe304 (magnetite) are also 
present. We have compared laboratory spectra of Fe203 
(Steyer 1974) to our observations. This oxide has strong 
bands in the 15-45 /tm region, that are sensitive to grain 
shape effects. We were unable to find reasonable fits to 
the observed spectra for spherical and for a CDE mixture. 
We conclude that Fe203 is not very abundant or absent 
in AFGL 4106. 
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Fe304 (magnetite) is quite featureless, except for a 
strong resonance at 17.5 /im and a weak one at 26.5 /im. 
These features are temperature dependent, not only in 
strength but also in shape and position and they have the 
tendency to shift to shorter wavelengths when the temper- 
ature decreases. At 150 K the 26.5 fim feature is split in 
two features at 25.4 and 26.7 /im, while the 17.5 /im fea- 
ture shifts to 17.1 /zm. Since we do not see strong features 
in the wavelength range from 17.0 to 17.5 /xm hardly any 
magnetite seems present. 

4.2.7. Unidentified features 

We find a narrow emission band at 47.8 /im; although the 
wavelength of this feature is close to a spectral structure in 
the responsivity curve of LWS and a instrumental artifact 
can therefore not be excluded, quite convincing evidence 
for its presence was given by Barlow et al. (1998) in the 
spectrum of several oxygen-rich sources with crystalline 
silicates at shorter wavelengths. According to the opti- 
cal constants of Steyer (1974), almandine (Fes AI2 (8104)3) 
peaks at 48 and 60 /im. However, other peaks of alman- 
dine, e.g. at 42 and 31.5 /im, are not observed. Therefore 
we exclude this possibility. We have not found a convinc- 
ing identification for this band, but given its width it is 
likely caused by a crystalline silicate. 

Several other peaks remain, which are not yet firmly 
identified (e.g. 16.8, 20.6, 22.9 and 24.1 fim). It is therefore 
(very) likely that more dust components are present. 

5. Radiative transfer modelling 

In this section we model the observed SED of 
AFGL 4106 with a one-dimensional radiative transfer code 
(MODUST, Bouwman & Waters, 1998; Bouwman et al. in 
prep.; de Koter et al., in prep.). From the structure seen 
in the N-band images it is clear that this approach is not 
entirely correct. Since AFGL 4106 is only marginally ex- 
tended for ISO, we assume that diffraction losses can be 
neglected. 

5.1. Dust species 

We have tried to collect optical constants of the dust 
species, that are identified in the previous section. Un- 
fortunately, in many cases the optical constants are given 
only for a limited wavelength range, and we had to com- 
bine data sets from different authors or extrapolate to 
cover the full wavelength range required for our modelling. 
The data sets used for the modelling can be found in Ta- 
ble |. 

We assume spherical grains and a power-law size dis- 
tribution and we use Mie theory to calculate the optical 
properties of the grain population. This approach allows 
us to calculate the emission for arbitrarily large grains. As 



we will show below, large grains are abundant in the shell 
of AFGL 4106. 

Since we identified both crystalline olivines as well as 
crystalline pyroxenes in our spectrum, it is likely that also 
both components are present in the amorphous state. For 
the amorphous olivines we applied the optical constants 
derived by Dorschner et al. (1995). Unfortunately they 
only published optical constants for two different Fe/Mg 
ratios (Mgo.8Fei.2Si04 and Mgi.oFei.oSi04). So, although 
the more Fe-rich olivine fits the data slightly better, it 
only gives a rough estimate of the Fe over Mg ratio for 
the amorphous olivines. 

Dorschner et al. (1995) also derived optical con- 
stants for amorphous pyroxenes ranging from MgSiOa to 
Mgo.4Feo.6Si03. Amorphous pyroxenes with the same per- 
centage of Fe as the amorphous olivines produce a peak 
around 21 /im, which is not seen in our ISO-spectrum. 
Only when the Fe content drops below 0.2 (i.e. Mg 
higher than 0.8), the peak vanishes. For Mg-rich amor- 
phous pyroxenes with a Mg content above 0.8 only two 
sets of optical constants are available (x=l and 0.95 in 
Mga;Fe(i_^)Si03). Since the emissivity of the Fe-poor py- 
roxenes is very dependent on the Fe content, it is difficult 
to give a reliable estimate of the contribution of the amor- 
phous Fe-poor pyroxenes based on these 2 sets only. Amor- 
phous pyroxenes, if present at all, will only marginally 
contribute to the IR-excess, due to their lower Fe content 
compared with the amorphous olivines, unless they have 
an efficient thermal coupling with Fe-containing dust, due 
to an aggregate or a core mantle structure. Since we have 
already found indications for separate grain populations, 
see e.g. the temperature differences found for the amor- 
phous and crystalline olivines, we assume that these two 
component grains do not exist and only the amorphous 
Fe-rich olivines are present. 

Crystalline forsterite has 3 different crystallographic 
axes. We found two studies of the optical constants of 
the 3 axes, one by Servoin & Piriou (1973) who used a 
synthetic sample and one by Steyer (1974), who used a 
natural sample. Both data sets do not cover the whole 
wavelength range. Although both sets of optical constants 
should describe the same material there were quite some 
differences. The Steyer (1974) optical constants were not 
able to reproduce the correct wavelengths of the 23.5 and 
33.5 /zm bands, found in laboratory transmission measure- 
ments (Koike et al. 1993; Jager et al. 1998), even when 
shape effects or core mantle structures were taken into 
account. The Servoin & Piriou (1973) data set predicts 
for spherical grains peaks at shorter wavelengths than ob- 
served. A CDE mixture produces peaks that are at slightly 
longer wavelengths than observed. It is likely that a some- 
what modified CDE mixture will give a good match. The 
difference between the Steyer and Servoin & Piriou data 
may be related to the origin of the samples, i.e. respec- 
tively natural and synthetic. Since Servoin and Piriou only 
have data longwards of 5 /im we extended the data set to 
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Dust species 


Chemical formula 




reference 


olivine 


Mgo.8Fei.2Si04 


0.1 - 0.2 


extrapolated 


(amorphous) 




0.2 - 500 


Dorschner et al. (1995) 






500 - 1500 


extrapolated 


forsterite 


Mg2Si04 


0.1 - 5.0 


Scott & Duley (1996) 


(crystalline) 




5.0 - 125 


Servoin & Piriou (1973) 


(all 3-axes) 




125 - 1500 


extrapolated 


water ice 


H2O 


0.1 - 1.25 


extrapolated 






1.25 - 333 


Bertie et al. (1969) 






333 - 1500 


extrapolated 


iron-oxide 


FeO 


0.1 - 0.2 


extrapolated 






0.2 - 500 


Henning et al. (1995) 






500 - 1500 


extrapolated 


corundum 


a-AbOa 


0.1 - 0.5 


extrapolated 






0.5 - 400 


Koike et al. (1995) 






400 - 1500 


extrapolated 


Table 5. The dust species used fo 


r the radiative transfer modelling 



the shorter wavelengths using amorphous forsterite data 
as pubhshed by Scott & Duley (1996), assuming that the 
absorption properties of the amorphous and crystalline 
grains do not differ much in this part of the spectrum. For 
a fair determination of the abundance of the crystalline 
forsterite we will apply Mie theory, since we use the spher- 
ical shape also for the other dust species. However for the 
best fit we apply CDE, since this fits the position of the 
peaks much better. In this case the abundance is decreased 
with 25% because of the influence of size and shape effects 
on the peak to continuum ratio of crystalline forsterite. 

To our knowledge no optical constants are avail- 
able for the three different crystallographic axes of the 
pure Mg-end member of the crystalline pyroxene solid 
solution serie, enstatite (MgSiOa). Jager et al. (1998) 
have measured ortho-enstatite with a small amount of 
Fe (Mgo.96Feo.o4Si03). Given the problems mentioned in 
Sect. 4.2.2 it is clear that we cannot reproduce all the 
peaks with the desired strength for this set of optical con- 
stants, e.g. the 32.8 iim peak. Therefore we decided not 
to include enstatite in our model fit. We will derive the 



enstatite abundance in a different way (see Sect. 5.2). 

An interesting point is how to combine the different 
grain components in model calculations. From an obser- 
vational point of view, there appear to be several separate 
grain components: the large temperature difference be- 
tween the amorphous and crystalline grains suggest that 
these components are not in thermal contact, which, as- 
suming that they are co-spatial, implies that they should 
be treated as two independent grain populations in the ra- 
diative transfer. The situation is less clear for the different 
crystalline silicates. They may very well be incorporated in 
the same grain, and this can change the optical properties 
of the material significantly. We have tried to produce a 
forsterite-enstatite-aggregate grain theoretically with the 
Bruggeman (1935) method, and although the results are 
promising we will not use these results since no good op- 



tical constants for the pure enstatite are available at the 
moment. For forsterite, with its three different crystal- 
lographic axes, we assumed an equal distribution of the 
orientation axes. This has been achieved by adopting 3 
different grain populations, one for each axis. The num- 
ber density, and thus the mass, of each population is one 
third of the total number density (mass) of the material. 

Our code is not well suited for large, non-spherical 
grains. Therefore for FeO we use instead of Mie a CDE 
shape distribution, which is not strictly correct for the size 
distribution we expect. Since shape effects also change the 
emissivity of FeO, the abundance of FeO derived from our 
model fit should be used with great caution. 

We have not included melilite, because optical con- 
stants with a sufficient large (0.1-100 /im) wavelength cov- 
erage were not available. 

The final point of caution is the fact that we used 
blackbodies to represent the stellar continua. This will in- 
fluence the result especially in the part where absorption 
from molecules dominates the stellar photosphere, but not 
for longer wavelengths where the dust emission dominates. 
A more detailed fit of the stellar photospheres will be pre- 
sented by van Winckel et al. (in prep.). 

5.2. The best fit model parameters 

Our best fit to the SED can be seen in Fi g. [l2| The model 
parameters of the shell are given in Table |6|, together with 
an estimate of their errors. Starting estimates for the inner 
and outer radius of the dust shell were derived from the 
TIMMI image, assuming a distance to AFGL 4106 of 3.3 
kpc. 



6.1. The dust 



This value will be discussed in Sect, 
shell is optically thin at most wavelengths. This implies 
that the inner and outer radius of the dust shell depend 
linearly on the distance, as does the duration of the mass- 
loss burst. The total mass in the shell scales roughly with 
the square of the distance, therefore the mass loss rate will 
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A [ycxm] 



Fig. 12. The SED of AFGL 4106 (thick hnc) and our best model fit (thin line) The input values are those from Table 0, 
with the exception that enstatite was not put in and for forsterite and FeO we applied CDE instead of Mie. 



also scale linearly while the density at the inner radius is 
inversely proportional to the distance. 

In Table |^ the dust parameters used in the model are 
given. To minimize the number of free parameters, we have 
assumed that all grain-species have the same size distri- 
bution. The typical errors in the grain size are 25%. In all 
models we have assumed a size distribution N{a) oc a~^'^ 
(Mathis et al. 1977). 



Apart from the uncertainty of the enstatite mass frac- 
tion, which influences the whole scheme, the errors in the 
mass fractions differ between the materials. For FeO the 
error is difficult to determine since the shape distribu- 
tion largely influences the outcome, an order of magni- 
tude cannot be excluded. The mass fraction of the amor- 
phous olivines is accurate to a few percent of the total 
mass. However the chemical Fe:Mg ratio is not well de- 
termined in the amorphous olivines, a change of 10% in 
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Stellar properties 



Lhot star 
Thot star 
Lcool star 
Tcool star 



1.3 X 10" 
7250 

7.4 X 10* 
3750 



K 

L© 
K 



Shell properties 



Inner radius 


4.2 


± 


1.0 X 10'" 


cm 


Outer radius 


4.1 


± 


1.0 X 10'^ 


cm 


p at the inner radius 


1.2 


± 


0.2 X 10"'* 


g/cm^ 


exp in p{R) = C X 7i°^P 


-2.2 


± 


0.15 




Mdust 


3.9 


± 


1.0 X 10"^ 


Mo 


<M> 


9. 


± 


2. X 10"* 


Mo/yr 



Table 6. The fixed stellar input parameters together with 
the Model parameters for the best fit, assuming a distance 
of 3.3 kpc. The errors for the shell parameters are internal 
model errors. 



Dust properties 



be interesting to compare this value with the olivine to 
pyroxene ratio seen for the amorphous silicates. However, 
because of reasons given in Sect. 5.1 this is not possible 
at the moment. 

Since we were not able to reproduce (in a single dust 
shell) the low temperature of crystalline H2O ice, the 
abundance of crystalline H2O ice is only based on the 
strength of the 43^m peak and is in that sense accurate 
to about 1% of the total mass. However, if there is a very 
cool H2O ice dust shell surrounding this object, which can 
be responsible for the LWS excess and missing strength of 
the 60 /im feature, the H2O ice abundance would probably 
increase at least an order of magnitude. 

We tried different models and only with large grains we 
were able to fit the shape, strength and position of the 10 
fim feature. Smaller grain sizes (spherical or not) would 
produce a distinct and relatively sharp peak at shorter 
wavelengths in contrast to the flat spectrum we observed. 



Grain size (/im) 

Mass fraction of crystalline forsterite (Mg2Si04) 
Mass fraction of crystalline enstatite (MgSiOa) 
Mass fraction of crystalline water ice (H2O) 
Mass fraction of crystalline corundum (Qf-Al203) 
Mass fraction of amorphous olivine (Mgo.8Fei.2Si04) 
Mass fraction of ironoxide (FeO) 

Table 7. The dust properties as derived from the radiative 
transfer model fitting. The ranges in the mass fractions 
come from the uncertainty about the enstatite to forsterite 
ratio. 



this ratio will only marginally influence the spectrum. The 
crystalline forsterite abundance is mainly based on the 
strength of the 23.6 and 33.6 micron peak and has an 
error of about 1% of the total mass. The abundance of 
enstatite is hard to determine since no optical constants 
were available for this material. So, we used an alternative 
way to determine the enstatite to forsterite ratio. From 
transmission measurements of enstatite and forsterite one 
may determine the wavelength dependent extinction coef- 
ficients. The 33.6 fim peak is caused by forsterite and the 
32.8 fim peak by ortho-enstatite. So the peak-strength ra- 
tio of these two features is indicative for the abundance 
ratio. Since the peak wavelengths of these two features 
are very close, temperature effects do not influence the 
results. If we use the transmission spectra published by 
Koike & Shibai, we find a ratio of enstatite to forsterite 
of 3. However, transmission spectra of forsterite published 
by Jager et al. (1998) give a peak strength of the 33.6 fim 
feature which is 3 times smaller than that of Koike & 
Shibai (1998), from which a ratio of enstatite to forsterite 
of 1 follows. We conclude that, given the large spread in 
laboratory data published so far, we find an abundance 
ratio of enstatite to forsterite between 1 and 3. It would 



0.4 - 6.0L he presence of these large grains causes a grey extinction 
4% at optical wavelengths. This influence s ou r previous de- 
4-11% termination of the extinction (see Sect. 3T) and therefore 
5% the luminosity of both stars. The luminosity of the stars 
17-15%y^,a,s derived using the standard interstellar extinction law 
70-65%^p]^^j,g et al. 1994) for both the interstellar and circum- 
"'^'^-stellar extinction. Since the colour excess changes much 
less when dealing with large grains, the stars are brighter 
than derived earlier in this paper. Using an iterative fit- 
ting scheme we re-derived the luminosity and interstellar 
extinction. We find a luminosity of 1.2 x IG^^d^L© for the 
hot star and 6.8 x 10^ (PLq for the cool star and an in- 
terstellar extinction of E{B — V) — 1.1 mag. Note that 
this value for the interstellar extinction is very close to 
the value which has been derived from the DIBs. Further- 
more, our model predicts an optical depth to the star of 
T = 0.6 at 6583 A, which is in very good agreement with 
the value t = 0.6 ±0.1 found by van Loon et al. (in prep.) 
for the optical depth at [N II]-line. Although r at 6583 
A is 0.6, the circumstellar E(B — V) is only 0.015 mag., 
due to the grey extinction (at optical wavelengths) of the 
large grains. So, the problem between the total extinction 
on one side and the circumstellar and interstellar dust ex- 
tinction on the other side is resolved by the presence of 
large grains in the circumstellar environment leading to a 
non-standard extinction law for the dust shell. 

Even with these large grains we were not able to fit 
the entire shape of the 10 ^m silicate feature. Especially 
around 13 fim our models predicted lower fiux levels than 
observed. When we subtract our model fit from the obser- 
vations, a 13 /im feature is evident. This feature resembles 
the 13 /Ltm feature observed in AGB stars with low mass 
loss rates, and is often attributed to corundum (q;-A1203) 
(Onaka et al. 1989; Begemann et al. 1997; Kozasa & So- 
gawa 1998). Therefore we also included corundum in our 
modelling, which indeed gives a much better fit around 13 
fim. 
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The total mass in the sheh, assuming a gas-to-dust 
ratio of 100 and a distance of 3.3 kpc, is 3.9 Mq. Taking all 
the different errors and assumptions (such as the spherical 
shape of the dust grains) into account, this mass may be 
overestimated by up to a factor two. 

While about 70 % of the dust shell consists of Fe- 
bearing amorphous olivines, about 7 to 15 %, depending 
on the assumed enstatite to forsterite ratio, constitutes of 
Mg-rich crystalline silicates. For the crystalline silicates 
only the end members of the olivine and pyroxene series, 
forsterite and enstatite respectively, are present. We note 
that the crystalline silicates only become important at 
the longer wavelengths. At 10 /im the amorphous silicates 
show a prominent peak but in this wavelength region there 
is no sign of the presence of the crystalline silicates. We 
conclude that the temperature of the crystalline silicates 
is lower than that of the amorphous ones. This temper- 
ature difference is confirmed by our model calculations. 
The temperatures at the inner radius are roughly 160, 125 
and 100 K, for respectively the amorphous olivines, crys- 
talline forsterite and the crystalline water-ice. These tem- 



peratures are higher than estimated in Sect. 4.2, but one 



should take into account that there is a temperature gra- 
dient and that the values mentioned here are maxima. The 
temperature difference is caused by the increased opacity 
in the UV, optical and near-IR when Fc is incorporated 
in the silicate. This implies that the two species are not 
thermally coupled, suggesting separate grain populations. 



5.3. Discussion of the model fit 

For some parts of the spectrum we did not obtain a satis- 
factory fit. The most obvious part is the LWS spectrum. 
We were not able to fit simultaneously the slope in the 
LWS part and the peak in the dust emission in the SWS 
part. A shallower density gradient would result in a broad- 
ening of the synthetic spectrum around 30 /im. We decided 
to fit the peak instead of the long wavelength part, also 
because of the differences between the 100 /im IRAS and 
LWS flux in this region, which required the steep density 
gradient. Note that the model spectrum convolved with 
the IRAS 100 /im sensitivity gives a flux of 125 Jy, which is 
lower than both the IRAS and LWS flux. The fact that the 
flux in our model is too low and that the LWS spectrum 
contains a prominent and very cool H2O ice band, gives an 
indication that an extra cool dust component, with likely a 
large amount of crystalline H20-ice, surrounding this ob- 
ject is required in the LWS beam. The presence of a large 
amount of crystalline H2O, which is generally not found 
in the ISM, suggests that this extra contribution is due 
to a previous mass loss phase. The discrepancy between 
IRAS and ISO, however, remains unsolved. GKH had a 3 
sigma detection of AFGL 4106 at 1.3 mm of 27 ± 9 mJy. 
In our model we find a flux of 17 niJy; this is within the 
errorbars. The ratio of 1.3 mm to (IRAS) 100 /im flux 



densities is reproduced by our model, indicating that the 
largest grains are significantly smaller than 100 /tm. 

Differences between the model and the broadband pho- 
tometry in the optical and near-infrared, where we see the 
photospheres of both stars, can probably all be explained 
by the fact that we used blackbodies instead of real stel- 
lar photospheres. A third discrepancy is between 4.5 and 8 
/im. We were not able to fit this flat part of the spectrum. 
Part of it can be explained by photospheric CO absorp- 
tion (around 4.5 fim). A few experiments with FeO seem 
to indicate that a proper treatment of the shape effects 
of FeO might produce a flat spectrum in this area. On 
the other hand, the TIMMI N-band image shows a bright 
central emission, which cannot only be explained by the 
photospheres of the stars. It seems that an extra (hot) 
dust component is present at the center and that its in- 
frared excess extends to the shorter wavelengths. Whether 
this represents a recent mass loss phase or a more stable 
circumstellar dust conflguration (circumbinary disk) is un- 
known. 

Previous modelling of AFGL 4106 has been performed 
by Volk and Kwok (1989) and HKV, who presented the 
same model, and GKH. Both models use in the center a 
single blackbody, with a temperature of 5000 (HKV) and 
4750 K (GKH), surrounded by a single dust component 
circumstellar envelope. Because of these big differences be- 
tween their and our model it is not useful to compare the 
inner radius and density at the inner radius found by them 
and us. 

It is interesting that HKV and the present study both 
flnd a density distribution which differs from p{R) oc i?^^' , 
with R the distance of the position in the dust shell to 
the central object. HKV found a steeper density gradient 
than we did, but this is within the error margins. The 
main difference between our parameters and the ones of 
GKH is the exponent in the density gradient. Since we 
have now spectral coverage longwards of 25 /im, we can 
exclude the possibility of a density gradient of p R^^. 
Their flt predicts a much broader IR excess around 25 
/im than observed. It should be noted that GKH fit the 
IRAS 60 and 100 /im points better than the present study. 
However, because of the significant broadening of the top 
of the energy distribution when applying a flatter density 
distribution, which is not observed, the IR-excess problem 
is more likely solved with an extra dust shell. The density 
gradient of p{R) ~ can be caused by two reasons: 

(I) The mass loss rate increased or (II) the outflow veloc- 
ity decreased with time. In both cases the density gradient 
increases due to the geometric expansion of the dust shell 
with time. So whatever the reason is for the present steep 
gradient, it is intrinsic but has steepened during the ex- 
pansion of the dust shell. 

Assuming different grain size distribution for the dif- 
ferent dust species might lead to an improvement of the fit 
and will result different mass ratios for the dust species. 
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A smaller size distribution for one dust component will 
result in a decrease of its mass ratio and vice versa. 

6. Discussion 

6.1. The distance and luminosity 

In the literature distances ranging from 0.67 - 3.3 kpc 
(Garcia-Lario et al. 1994) are given. However one should 
note that AFGL 4106 is located in the direction of Ca- 
rina, so we are looking along a spiral arm, which makes it 
hard to relate column density and distance. Fortunately, 
we have spatial information from the N-band image. The 
inner radius of the dust shell is about 1.", and together 
with the expansion velocity found in CO (Josselin et al. 
1998) and [Nll]-line (van Loon et al., in prep.) of 30 to 40 
km/s and the present-day shape of the SED, one can put 
a lower limit to its distance. Assuming that the outflow 
is spherical, which seems plausible according to our mod- 
elling results and the [N II]-line velocity distribution by 
van Loon et al. (in prep.), and that the velocity has not 
changed significantly in time, one can derive a kinematical 
age for the dust shell of 135 years times the distance in kpc. 
Due to the expansion of the dust shell it is expected that 
the circumstellar dust extinction becomes smaller, lead- 
ing to an increase in the brightness at short wavelengths. 
AFGL 4106 has already been identified in the Cape Obser- 
vatory Photographic Durchmusterung from 1895 to 1900 
(Gill and Kaptcyn 1900) with a photographic magnitude 
Mpg = 9.9 mag. In the Cordoba Durchmusterung (Thome 
1914) a visual magnitude of 9.4 has been estimated for 
AFGL 4106. In the extended Henry Draper Catalogue, 
with observations taken between 1922 and 1937 (Can- 
non and Mayall 1949), its photographic magnitude was 
10.0 ± 0.3. On 18 March 1988 Hrivnak et al. (1989) mea- 
sured a B magnitude of 10.24 ± 0.01 (Mpg = Mb - 0.11) 
and a V magnitude of 8.73±0.01. From these observations 
one can derive that AMpg = 0.2 ± 0.5 mag. over the last 
90 years and AMy = —0.7 ± 0.5 mag. over the last 75 
years, assuming an error of 0.5 magnitude in both Durch- 
musterungen. So the visual and photographic band trends 
contradict each other, although the errorbars are signif- 
icant. Since the hot star is the dominant component at 
these two wavelengths it is possible that the star lowered 
its temperature, from a temperature with its peak lumi- 
nosity at the photographic band to a temperature with 
its peak luminosity at the visual band. However, this is in 
contradiction with the present temperature, which peaks 
around the photographic band and is also not expected 
from evolutionary considerations. 

In any case, the stability of the optical brightness of 
AFGL 4106 over a timescale of w 90 years can be used 
to constrain its distance and luminosity. Using our dust 
model fit we can calculate the expected SED of the system 
and the visual and photographic brightness in the past, 
assuming that the dust composition, the outflow velocity, 
the temperature and the radii of the stars did not change. 



For the photographic magnitude, it is safe to assume that 
the star has not brightened more than 0.3 magnitude over 
the last 90 years. This puts the star at a distance of at 
least 3.15 kpc and probably further away. On the other 
hand the brightning in the V-band over the last 75 years 
of AMy = — 0.7±0.5 puts the star at a distance of 1.51q;4 
kpc. There is only a very small overlap region between 
these two trends, suggesting a distance of 3.3 kpc, just 
at the upper limit quoted by Garcia-Lario et al. (1994). 
The size of this overlap region is mainly based on the 
estimated errors of the Cape Observatory Photographic 
Durchmusterung and the Cordoba Durchmusterung. 

This distance of 3.3 kpc implies that the luminosity 
is 1.3 X 10^ and 7.4 x 10'' Lq for the hot and cool star 
respectively, again pointing to a massive binary system. 
Their luminosities point to stars with a main sequence 
mass between 15 and 20 Mq, using evolutionary tracks by 
Maeder & Meynet (1988). Other properties of AFGL 4106 
also support the high masses of both stars. The expansion 
velocity of 30 to 40 km/sec is unusually high for AGB 
stars, and the relative large grains are typically found in 
the outflows massive supergiants (Jura 1996). 

It is likely that the warm star is responsible for the 
expelled dust shell, and is now in the very rare evolution- 
ary phase of the post-Red-Supergiants, which makes it a 
blue supergiant or WR star progenitor. Only a few other 
stars are in the same evolutionary status, e.g. IRC-I- 10420 
(Jones et al. 1993), and the central star of the radio nebula 
G79.29-h0.46 (Trams et al. 1999). 

6.2. The dust shell 

At a distance of 3.3 kpc the mass loss episode took about 
4.3 X 10^ years and stopped 450 years ago. During this 
period the star expelled 3.9 Mq, assuming a gas to dust 
ratio of 100. This implies that the average mass loss rate 
is roughly 9. x 10"'* Mq/jt. 

In our modelling we assumed spherical symmetry, but 

we do see substructure in our N-band images. This is also 
expected, since it is likely that the M-type supergiant will 
influence the spatial distribution of the dust and gas, ex- 
pelled by the hot star, by transferring orbital momentum 
to the gas and dust. This effect can create a density en- 
hancement in the plane of the binary. This possible en- 
hancement can lead to an over- or under-estimate of the 
total mass, depending on the angle at which we observe 
this system. However, the shape of the [N Il]-line derived 
from different rotation angles of the slit, suggest that the 
shape of the dust shell is close to a sphere (van Loon et 
al., in prep.). This implies that the orbital seperation be- 
tween the two stars is probably relatively large and there- 
fore the influence of the companion on the mass-loss not 
so important. So, although the spherical symmetry and 
other assumptions used in our model calculations might 
be somewhat simple, it is not expected that the mass-loss 
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rate and the total mass found in our calculations would 
change more than a factor 2. 

6.3. Formation of crystalline silicates 

The formation of crystalline silicates is still largely un- 
known. Several researchers have tried to quantify the con- 
densation sequence of oxygen rich dust (e.g. Tielens, 1990; 
Gail, 1998). However, these are mainly based on the differ- 
ent condensation temperatures, therefore assuming ther- 
mal equilibrium; it is not known if this assumption is valid. 
Condensation of dust species and chemical reactions will 
only occur if: (I) the density is high enough and (II) the 
temperature is suitable (low enough for condensation and 
high enough for chemical reactions). Since both the den- 
sity and the temperature decrease with the distance from 
the star, it is likely that at a certain moment the dust 
structure freezes out. It is this stable dust configuration 
which we sec at the moment. Taking the condensation 
sequences mentioned above, it is expected that the first 
silicate that will form is forsterite. When the temperature 
becomes lower and the density is still high enough this 
will transform into enstatite. Both species are expected 
to form above the glass temperature and are therefore 
thought to be crystalline. We see evidence in our spec- 
trum for both forsterite and enstatite. This implies that 
the forsterite to enstatite transition is stopped before the 
forsterite was completely transformed. There might be 3 
reasons for this: (I) The density became too low, (II) the 
temperature dropped too rapidly or (III) other chemical 
reactions took over. It is likely that it is a combination of 
these three. One chemical reaction that could take place at 
lower temperatures is the incorporation of Fe into the sili- 
cates. Since our spectra show that the crystalline silicates 
are very Mg- rich and the amorphous silicates contain a 
lot of Fe, it is likely that this incorporation of Fe results in 
a destruction of the crystal structure. The clear chemical 
separation between the amorphous and crystalline mate- 
rials is intriguing. It appears that the inclusion of Fe in 
Mg-rich silicates is a kind of runaway process. 

Tielens et al. (1998) proposed a scenario for such a 
chemical separation between amorphous and crystalline 
silicates. When the temperature becomes low enough Fe 
may react with the Mg-rich crystalline silicates. The opac- 
ity will increase, due to the incorporation of Fe, and there- 
fore the temperature of the grain. This process will act as 
a thermostat, incorporating just sufficient Fe in the grains 
to keep the temperature near 800 K where Fe can just 
diffuse in. This temperature is below the glass tempera- 
ture for Fe-bearing silicates and the lattice cannot attain 
its energetically most favorable structure, thus leading to 
an amorphous structure. Because of the increased opacity 
due to the incorporation of Fe, it is likely that the temper- 
ature of the other grains behind this Fe-reaction zone will 
decrease, leading to an even higher difference between the 
Fe-rich (amorphous) and Mg-rich (crystalline) silicates if 



thermal coupling between gas and dust can be neglected. 
Grains that already contain a small amount of Fe will 
be hotter and therefore the reaction rate for the adsorp- 
tion of Fe will be higher then for grains without any Fe. 
This would lead to the required runaway process. If this 
scenario is correct, it would imply that the crystalline sil- 
icates are the primordial condensation products, and this 
will give us new insights in the dust nucleation. Although 
this would lead to the required chemical separation be- 
tween the crystalline and amorphous materials, it is not 
clear why the amorphous olivines should contain more Fe 
than the amorphous pyroxenes, if present at all. 

Besides the forsterite-enstatite condensation sequence, 
Tielens (1990) also presented a condensation sequence 
starting with corundum. The next condensation product 
would then be melilite, for which we find some evidence. 
According to the same scheme diopside would form. Since 
certain strong features of this material are not found in 
our ISO-spectrum, this may imply that the dust forming 
process along this condensation sequence line froze out 
around amorphous melilite. 

Gail and Sedlmayr (submitted to A&A) applied non- 
equilibrium calculation for the dust formation in outfiows 
of M stars. They do explain the formation of several dust 
species, such as the amorphous Fe-rich olivines, however 
they do not consider, based on laboratorium experiments, 
the formation of pyroxenes and are not capable to ex- 
plain the co-existence of pure forserite particles and Fe- 
rich amorphous olivines. 

7. Conclusions 

The main results of this study can be summarized as fol- 
lows: 

1 AFGL 4106 is a (double-lined spectroscopic) binary, 

consisting of a warm star with Toff = 7250 K and a 
cool companion with T^g = 3750 K. 

2 The luminosity ratio of both stars Lwarm/Lcooi = 1-8, 
which indicates that both stars are evolved. The warm 
star is likely responsible for the dust shell. 

3 The expansion velocity of the shell, the grain size dis- 
tribution and the minimal variations in the photo- 
graphic and V band over the last 90 years all suggest 
a high luminosity. We conclude that the stars are at 
3.3 ± 1.0 kpc and have L/Lq = 1.3 x 10^ and 7.4 x lO" 
respectively. These values correspond to main sequence 
masses of 15 to 20 Mq. The warm star is probably 
evolving rapidly to the blue part of the HR-diagram 
and may evolve to a blue-supergiant (ending its life 
like SN 1987A) or to a Wolf-Rayet phase. 

4 The main component of the dust are large (« 1/im) 
amorphous, Fe-rich olivines. If also amorphous pyrox- 
enes are present, they will have a much lower Fe con- 
tent. We find a rich spectrum of narrow solid state 
emission bands in the ISO-SWS and LWS spectra. 
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which we identify with crystaUine ohvines and pyrox- 
enes. These grains are Mg-rich and Fe-poor and have 
an abundance between 7 and 15 % by mass, depend- 
ing on the assumed cnstatitc to forstcritc ratio. We also 
find evidence for the presence of FeO, AI2O3, meUUte 
and crystaUine H20-ice. Shape effects can have an im- 
portant influence on the derived abundances. 

5 The temperatures for the different dust species are 
quite different, this is caused by variations in their 
absorption characteristics at visual and near-IR wave- 
lengths. These temperature differences imply that the 
different dust species are not thermally coupled, di- 
rectly (e.g. as a composite grain) nor indirectly (e.g. 
via gas-dust interactions). 

6 The crystalline silicate bands can be divided into nar- 
row (FWHM/A < 0.03) and broad (FWHM/A > 0.03) 
features. The wavelengths of the broad features match 
well with forsterite, while those of most of the narrow 
features line up well with enstatite. The width of the 
features is significantly smaller than those observed in 
laboratory spectra. The abundance ratio of enstatite 
to forsterite is between 3 and 1. 

7 The mean mass-loss rate was « 9. x 10^'' M0/yr for a 
period of about 4.3 x 10^ years and probably increased 
during this period. The mass loss stopped 450 years 
ago. The total mass expelled, assuming a gas to dust 
ratio of 100, is 3.9 M©. 

8 We find a cool dust component, which cannot be fitted 
by our dust model. This component is either an older 
mass loss phase, or an incorrect background subtrac- 
tion of the ISM. We also found indications from the 
spectrum between 4 and 7 /xm and the 10 /xm imag- 
ing that a third dust component, close to the binary 
is present. Whether this represents a recent mass loss 
phase or a more stable circumstellar dust configuration 
is not known. 

Because of its (IR-)brightness AFGL 4106 will be a key 

object in the further study of (crystalline) dust forma- 
tion. High resolution imaging and spectroscopy both in 
the optical and IR will allow us to constrain the circum- 
stellar versus interstellar extinction, the photospheres of 
the components, density and abundance gradients through 
the CSE, and thus the conditions for dust formation and 
evolution. 
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